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Abstract 

Symmetrization of bosonic wave functions produces bunching and low temperature 
phenomena like Bose-Einstein condensation, superfludity, superconductivity, as well as 
laser. Since probability is conserved, one might expect such bunchings at one place 
would lead to a depletion or a destructive interference at another place. This indeed 
takes place in high-energy scatterings, leading to cancellations between Feynman dia- 
grams with permuted (real or virtual) boson lines. It is such cancellation that allow 
the Froissart bound to be obeyed, and the meson-baryon amplitudes in large- N c QCD 
to be consistent. These lectures review the derivation and some of the applications 
of nonabelian cut diagrams, which are resummation of Feynman diagrams, but with 
these interference effects explicitly incorporated. Unlike the traditional method where 
cancellations occur between Feynman diagrams, making it necessary to compute each 
diagram to subleading orders, destructive interference has been taken care of from 
the start in nonabelian cut diagrams. No further cancellation takes place, so each 
nonabelian cut diagram may be computed only in the leading order. 
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1 PRELIMINARIES 



1.1 Introduction 

The symmetrization of bosonic wave functions leads to an effective attraction. At low tem- 
peratures this constructive interference gives rise to the Bose-Einstein condensation, super- 
fluidity of Helium, and superconductivity, as well as laser. 

Probability is conserved. If there is a concentration of wave function at one place then 
there is presumably a depletion at another. How can one observe this depletion, or destruc- 
tive interference? It turns out that this can be seen in high energy processes. The effect is 
not as dramatic as the low-temperature constructive interference processes, but its presence 
is absolutely necessary to render theoretical consistency, as we shall see. 

A general discussion of this interference can be found in Sees. 3 and 4. The interference 
effect can be incoporated in our calculations by using the the multiple commutator formula [|lj 
and the nonabelian cut diagrams |2j . These cut diagrams may be considered as resummations 
of Feynman diagrams where the interference of the identical bosons in the diagrams, whether 
real or virtual, has been explicitly accounted for at the very first step, the step of drawing the 
diagrams before any computation begins. With the interference thus included, subsequant 
cancellations are not needed and will not occur. In contrast, if as usual Feynman diagrams 
are used, substantial amount of delicate cancellations between diagrams will occur at a 
later stage as a result of such destructive interference. This is often much more difficult to 
calculate as will be elaborated below. 

A fairly dramatic example of this can be found in the QCD theory with a large number 
of color (N c ^> 1). Its n-meson tree amplitude in the one-baryon sector, viz., meson-baryon 

scattering producing n— 1 mesons in the final state, depends on N c as Nc 2n @- In particular, 
taking n = 1, this shows that the meson-baryon Yukawa coupling constant go like \/W c . 
Since the baryon propagators are of order unity, individual Feynman tree diagrams will go 

like Nc , which is n — 1 powers of N c larger than what the full amplitude should be. Unless 
a huge cancellation occurs when we sum up the n\ Feynman diagrams, the theory will not 
be consistent @]! Such cancellation indeed occurs, it is a manifestation of the destructive 
interference mentioned earlier, but to see how that works directly using Feynman diagram 
is very difficult task. Not so with the nonabelian cut diagrams || because the interference 
effects leading to the cancellations have already been built in. We will not pursue this subject 
further here because it is a bit outside of the scope of this workshop. 

What we will discuss is destructive interference in high-energy elastic scattering taking 
place in loops, where identical bosons are present in the intermediate states for the inter- 
ference to take place. Individual Feynman diagrams with many loops often contain high 
powers of f m 2//^du/u) ~ Ins [§, [7|, @. If not cancelled, the Froissart bound would be 
violated. When Feynman diagrams are summed, such offending In s do get cancelled |J 
as a result of the destructive interference mentioned earlier. However, it is difficult to see 
how this happens in traditional calculations for two reasons. First, one has to do a lot of 
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virtual work to compute terms that eventually get cancelled. Second and more seriously, 
such computations can usually be carried out only in the leading-log approximation. If these 
leading powers of Ins are cancelled out, then one needs to compute the sub leading, the 
sub-subleading (etc.) terms in order to obtain a non-zero finite result at the end. Such tasks 
are often impossible except in very low orders. Again, the nonabelian cut diagrams bypass 
these difficulties [§, [TO], |llj because the destructive interference effect has been taken care of 



before any computation even begins. 
1.2 Conventions 

We shall adopt the following conventions for the metric g^ u , the Dirac matrix 7 M , and the 
Dirac spinor u: 

<T = (+ )diag 

{7^,7^} = 2<r 

UU = 2M . (1.1) 

The plane-wave states are normalized covariantly, 

(p'\p) = (2n) 3 (2p°)5 3 (p'-p) , (1.2) 

the S'-matrix is related to the covariant T-matrix by the formula 

% = l /4 + (2 7 r)^ 4 (---)T / ,, (1.3) 

where <5 4 (- • ■) is the overall energy-momentum conservation 5-function, and the unitarity 
relation is given by 

- % - (T fi - 7>* ) = Im (T fi ) = ~ £ / T* fn d® n T ni , (1.4) 
in which the expression for the n-body invariant phase space is 

For elastic scattering with initial momenta pi and final momenta p' { , momentum conserva- 
tion dictates that P1+P2 = p[+P2- The Mandelstam variables are s = (P1+P2) 2 , t = (pi—p[) 2 , 
and u = (pi — p'2) 2 ', they are related by s + t + u = 2(Mf + M|). y/s is the centre-of-mass 
energy, and A = \f—i is the momentum transfer. We shall be concerned with situations 
where s is much larger than —t and Mf. The elastic cross-section is related to the elastic 
amplitude by 

dff = /l| T/ ,i^ 2 = _^/ WA , (1 . 6) 

and the total cross-section is given by unitarity to be 

atot = -Im(Tii) . (1.7) 
s 
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1.3 Feynman Rules 

The Feynman rules for the T-matrix elements that will be needed are the following. 

• Quark propagator: (M + r )-p)/{p 2 — M 2 + it) 

• Gluon propagators in the Feynman gauge: —g^ u / (p 2 + ie) 

• Quark-gluon vertex: gt a ^ a 



• Triple-gluon vertex: g{if a bc){g af3 {Pi - P2) 1 + 9 M {Vi ~ PzY + g ia {Pz ~ Pi) 13 } 

• Four-gluon vertex: g 2 {f ab efcde{g ai g fiS - g a6 g Py ) + facef b de(g aP g' rS - g aS g Py ) + 



• External wave functions: e* x (p) , e\(p) , u\(p) , u\(p) 

• Loop integration: % j d i q/{2^) A 

• Other factors: an overall minus sign, a minus sign per fermion loop, and symmetry 



factors. 



1.4 Color Algebra 

The SU (N c ) color matrices t a obeys the commutation relation 




(1.8) 



The structure constants satisfy the sum rules 



fabcfabd ^c^cdi i fadgfbedfcge i „ fabc — ^cfabc • 



(1.9) 



See Fig. 1 for a graphical representation of (1.8) and (1.9). 
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(a) 



(b) 



a b b a 



a b 



a b 

c 



(c) 



ab ba ab ab 



(d) 



= 2c 



(e) 



Figure 1: Graphical representation of the color-matrix relations (1.8) and (1.9). 



2 HIGH-ENERGY ELASTIC SCATTERING 



2.1 s-Dependence and Regge Poles 



A(p;) 



B( P ;) 



A(p t ) 

<— s 
B(r) 



Figure 2: An elastic scattering process A + B ^ A + B. 



Consider the tree diagram in Fig. 2 for the elastic process A + B — > A + B. For simplicity 
we shall assume A and I? to have the same mass. If all the particles are scalars then the 
T-matrix amplitude is —g 2 /t and is energy independent. On the other hand, if the exchange 
particle C is a photon, then the vertices add on an additional factor ~ — (2pi)-(2p 2 ) — —2s, 
now the amplitude grows linearly with s. In general, the amplitude grows like s e if the spin 
of the exchanged particle C is i. To see that, analytically continue the amplitude to the 
physical region of the t-channel process B + B — > A + A, where t > 0, s < 0, u < 0. In its 
CM system, the total angular momentum is £, so the amplitude is proportional to Pe(cos6 t ), 
where 9 t is the CM scattering angle for this process and it is related to s and u by 

s = -^(l-cos0 t ) , u= -|(l + cos0 t ) . (2.1) 

We can now analytically continue back to the physical region of the s-channel process and 
examine the limit when s — > oo. Since 

P t (z)~J, (Re£>~) (2.2) 

for large z t = — cos t ~ s, the amplitude grows like s £ as claimed. 

What if we have a complicated diagram so that the t-channel angular momentum is not 
fixed? Then the amplitude is given by a partial-wave expansion 

T(s,t) = Y,a e (t)(2£ + l)P e (z t ) , (2.3) 

i 

say with some 'mean' angular momentum 1 = a(t). If the spread of angular momentum 
is small, one might expect simply by interpolation that the amplitude still grows like s a ^\ 
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L +i°o 




Figure 3: Various contours used for the Sommerfeld- Watson representation. 

This naive expectation will be justified below even without this narrow-spread assumption. 
The function a(t) is known as a Regge trajectory. 

Let us assume ae(t) to be an analytic function of £, dying off sufficiently fast at infinity 
to enable the Sommerfeld- Watson representation to be used. We can then replace the sum 
(2.3) by a contour integral wrapping around the positive real £-axis, 

r( s , f ) = -|/ r ^ + i)« (W ^, (2.4) 

as shown in Fig. 3. 

Suppose a e (t) has (Regge) poles located at £ = ««(£) with residue $(t) / (2ai(t) + 1). 
Opening up the contour T and moving it to the left of the poles, we get 

T(s,t) = -^ / L+ ^^+l)a^)^^-^-^^P Qi(t) (-^ , (2.5) 
1% Jl-ioo sm7rl: Y sin7ra!j(Ej 

where L is smaller than the real part of all ati(t). The high-energy behavior is now dominated 
by the rightmost Regge pole a(t) to be 

{-z t ) a{t) ~ (2s/t) a{t) , (2.6) 
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agreeing with the naive expectation. 

In trying to get to this quickly I have glossed over a number of fine points. I have 
not discussed why T(s, t) is analytic but that can be established at least in perturbation 
theory. I have ignored the spins of the colliding particles which complicate matter but pose 
no essential difficulty Helicity is conserved at high energies and the dominant amplitude 
is helicity independent, so for elastic scattering everything is essentially the same as if the 
beam and target particles were without spin anyway. Another caveat is that the asymptotic 
behavior (2.2) is valid only when Re(£) > — | so one might worry about the general validity 
of (2.6). It turns out that this can always be fixed up to render (2.6) valid. 

I have also not discussed why it is reasonable to assume ae(t) to be analytic in £. In fact, 
this is not even true unless the signature factor is taken into account. 

To see what signature is let us suppose the particle A is its own antiparticle: A = A. 
Then there is a forward-backward symmetry for the t channel process B + B ^ A + A so 
that the amplitude T(s,t) is s <-> u symmetric. From (2.1) and (2.3), we see that only even 
angular momenta £ are present; ai(t) would be zero for odd £. One would therefore expect 
that only even £ values are smoothly connected by the function ai(t). 

This example suggests that we should always decompose an amplitude T(s, t) into its 
s <-> u symmetric/antisymmetric parts T±(s,t). Both a±e(t) are analytic but they may have 
different singularities, so the Regge trajectories for these two will generally be different. 
These two kinds of Regges trajectories are said to have even/odd signatures. When we 
replace T in (2.3) by T±, the right-hand side should be replaced by Pe(z t ) ± Pe(—z t ), so the 
last factor in (2.5) and (2.6) become 

1 [p ai(t) (-z t )±p ai{t) (z t )} - * ,A (-ztr®±(ztr®] 



sin ~n,(/ ) sin ~o,(/) 



t) sin7ra!i(t) 



[1 ± e™ l{t) ] . (2.7) 



The factor (1 ± e™^)/ sin7ro;(t) blows up at even/odd integer angular momenta J. This 
pole at t — rrij corresponds to a resonance in the t-channel with mass mj and spin J, so 
the Regge trajectory a(t) may be regarded as a trajectory connecting particles in the J 
- mj plane. Extrapolating this trajectory to negative t, it will determine the high-energy 
behavior of the s-channel amplitude, and this is the power of the Regge theory! However, 
this would work only if the singularities of ae(t) in £ are indeed poles, an assumption that 
requires explicit verification and is not automatically guaranteed. In fact, as we shall see, 
the pole assumption appears to be true for the gluon trajectory but not the Pomeranchuk 
trajectory that is discussed see below. B Experimentally the total cross-section grows with 
energy, thus from (1.7) and (2.7) it follows that the leading regge intercept is at least 1, 
cc(0) > 1. However, the Froissart bound forbids the total cross-section to grow faster than 
(Ins) 2 asymptotically, so this leading singularity, called the Pomeranchuk trajectory or the 
Pomeron for short, must have an unit intercept a(0) = 1, but then it cannot be a simple 
pole for otherwise the total cross-section would be a constant. It may be a double pole or 
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something more complicated. The usual view is that in QCD it is a composite object made 
up of two and more reggeized gluons. 



2.2 QED to 0{g 6 ) 

Let us now turn to actual calculations of high energy processes in QED. We will follow the 
book of Cheng and Wu [6] from which one can find more details, as well as references to the 
original literature. We will concentrate specifically on electron-electron elastic scattering in 
the following. 

The incoming beams will be assumed to move in the z-direction. At high energy it is 
convenient to use the lightcone coordinates a± = a°±a 3 , in terms of which the dot product of 
two four-vectors a = (a+, a_, a±) and b = (b + , b±) is given by ab = j(a + b^+a^b + ) — a±-b±. 
The momenta of the incoming beams in the CM system are p\ = {\/s, 0, 0) and p 2 = 
(0, y/s, 0), with masses ignored. We will make it a practice to draw pi on top and p 2 at the 
bottom in all the scattering diagrams. 

We shall label the loop momenta by and define the scaled '— ' momenta to be Xi = 
qi-/\fs. The measure for loop integration in lightcone variables is 



(27T) 



rA = 



8tt 2 
1 

47T 



dq+dq ~ (2^ V = I 
[Dq + ]dx[Dq ± ] . 



-2ni 



dx 



(27T) 



;d 2 q± 



(2.8) 



The momenta along the upper electron line are Pi + Q, with Q being some combination 
of qi (see Figs. 4 and 5 for illustrations). In the leading-log approximation when |A| = v / — t 
is fixed and s — > oo, we can approximate the electron propagator by 



M + 7(pi + Q) 



(pi + Q) 2 - M 2 + ie 2 Pl -Q + ie 



(2.9) 



Since 



M + 7P1 = 5Z u ^(pi)«a(Pi) > «a'(pi)7 m «a(pi) = (2^)^'. 



(2.10) 



we may simply drop all the 7 matrices and spinor numerators, and consider the vertices to 
be 2 P i. The same is true for the lower electron line where vertices can be taken to be 2p 2 . 
Physically, this expresses the simple fact that the large electric current carried by an electron 
at high energy is a result of its fast motion, with its magnetic moment contributing only a 
negligible amount. For the same reason the current of a spin-1 particle can be taken to be 
2p M , which is valid even when that is the color current imbedded in a triple- gluon vertex in 
QCD, because the other two terms in the vertex are of order A/y^ and can be neglected. 



The procedure to follow to compute the high-energy amplitude is as follows: 
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1. Use residue calculus and flow diagrams to compute the '+' integration. 

2. Then obtain powers of (Ins) from the '— ' integrations. 

3. The qi± integrations are never explicitly carried out. 

2.2.1 The '+' integration 

Let us consider step 1 in more detail. 

First of all, a-b is linear in a+ and b+, hence the denominators of all propagators are 
linear in q i+ . This then provides one simple pole per propagator for the '+' integrations. 
The integration contour is originally along the real q i+ axis, but we will close it with an 
extra half circle on the upper plane or the lower plane. By using — q i+ as the integration 
variable if necessary, we can and will always choose to close the contour in the lower plane. 
The result of the integral is therefore equal to —2-ki times the residue at the pole, summed 
over all poles in the lower plane. The important question to decide then is which are the 
poles in the lower plane, and which are the poles in the upper plane. This clearly depends 
on the sign of the ' — ' components, but they are integration variables and hence vary. The 
flow diagrams are invented to get the various possibilities sorted out [6]. 



Pi 



R-A 



R-q 



3 I 



(a) 



(b) 



Figure 4: A one- loop Feynman diagram (a) and its flow diagram (b). 

Figs. 4(a) and 5(a) are Feynman diagrams and 4(b), 5(b), and 5(c) are flow diagrams. 
The latter are diagrams indicating flows of the '— ' components of the momenta, subject 
to momentum conservation at each vertex. Arrows going one way (clockwise or counter 
clockwise) correspond to positive '— ' components, and arrows going the other way correspond 
to negative ' — ' components. Which is which does not matter because we can always change 
the integration variable from q i+ to — q i+ to compensate our choice. 
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Pi+A A+q, Pj+q, P, 

•< | i -■s i -< 

i i i 

qr A i qr^i i q, 

i i i 

i i i 

i i i 

p-a" v 2 -% V 2 -<\ " P 2 
(a) 

7 6 It 



: : : x : : 

5) 4) 3 1 5) 4} 3 1 

i i X iii 



2 1 2 1 

Figure 5: A two- loop Feynman diagram (a) and its flow diagrams (b) and (c) 



Note that we have not included in Figs. 4 and 5 any flow diagram in which the arrows 
go around in the same direction in a close loop, for this gives rise to poles all in the same 
half plane so the '+' integration would vanish. In making up flow diagrams, such closed 
loops are always excluded. In particular, this means that at the injection vertex where P2 
enters, and at the ejection vertex where p' 2 leaves, the arrows are always arranged in the 
way shown in these diagrams. This rule determines uniquely the flow pattern of a one-loop 
diagram, as shown in 4(b), but it does not do so for diagrams of two or more loops, as the 
flow direction at the boundary of two loops can always be reversed. For two loops we have 
two flow diagrams, 5(b) and 5(c), and for more loops there are more flow diagrams. 

Since we would like to have as few poles as possible in the lower plane, we would choose 
those directions in every loop with less arrows to correspond to poles in the lower plane. 
These are indicated by a cross (x) in the flow diagrams. 

In the case of the right-hand loop in 5(b), there are two arrows along each direction, so 
it would seem to have two equally valid choices of poles for the lower plane. However, the 
ones indicated are preferred for two reasons. 

First and foremost, it is never a good idea to choose a pole at the boundary of two loops. 
A pole within loop i determines the value of q i+ by the location of the pole, but a pole at the 
boundary of loops i and j determines only the combination q i+ + qj + which is less convenient. 

There is a second reason for that choice. With the approximation (2.9), q i+ never appears 
on a propagator along the top line, so apparent poles there are actually absent. With that, 
the choice indicated actually has only one pole per loop. 
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2.2.2 The '-' integration 



The inverse propagators along the upper electron line is given by (2.9) to be sX, where 
X = Q_/y/s is some linear combination of Xi, and each Xi in a flow diagram is by definition 
positive. The Ins factors are typically obtained from J e dxi/xi, where e ~ A 2 /s represents a 
cutoff of Xi above which the approximation (2.9) is valid. For that reason, we can usually 
ignore all X{ compared to 1 in the leading-log approximation. 

To proceed further it is convenient to adopt some shorthand notations. We will write 
a m = ^ is the four-momentum of the mth line. If its pole is in the upper plane, we 
let o?" 1 to be its propagator. If it is in the lower plane, we let d" 1 to be its residue multiplied 
by y/s. It is understood that all of them are evaluated at the lower-plane pole positions, and 
the minus sign in the gluon propagator is included in this factor. 

From (2.8), we see that the T-matrix after the Dqi + integrations is given by the formula 



Here t is the number of loops in the diagram, D = ]J m d m , and N is the numerator factor 
coming from products of vertex factors. The lower limits of Xi can be taken to be e = A 2 /s. 
The minus sign comes from the overall sign factor in the last item of Sec. 1.3. 

2.2.3 Fig. 4 

From Figs. 4(a) and 4(b), we see that the pole is located at q + = a\j (l — x)y/s ~ ai/y/s, with 
d± — (1 — x) ~ 1. The rest of the (i-factors are d 2 = a 2 , d 3 = a 3 , c? 4 = sx, so D ~ l-a 2 -a 3 -sx. 
The numerator is iV = g 4 [(2p 1 )-(2p 2 )] 2 = (? 4 (4s 2 ), so the T-matrix of this diagram is obtained 
from (2.11) to be 




(2.11) 



T 




-(ln S )// 2 (A) , 



(2.12) 



where 




(2.13) 



In terms of convolution defined by 




(2.14) 



the function 



/„ is the nth-power convolution of I\ with itself, which we write as 



J B (A) = (♦AHA) . 



(2.15) 
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The impact-parameter-space function F(b) is the Fourier transform of F(q±), 

Fib) = J 0- 2 F(q ± )exp(-tq ± -b) . (2.16) 

Such functions are convenient because in this space convolutions turn into ordinary products 
in the usual way: [F*G](6) = F(b)G(b). Hence I n {b) = [h(b)] n . 



2.2.4 Fig. 5 

First consider flow diagram 5(b). The poles are located at q 2+ = a 2 /(l—x 2 )y/s ~ a 2 /-y/i, with 
d 2 — 1, and q i+ = a^/xi^/s with d 3 = x\. The rest of the d- factors are d\ = (p 2 — qi)+(p2 — 
qx)--ai ~ -a 3 /xi-ai ~ -a 3 /xi,d4 = -(gi-g 2 )+(gi-g2)-+04 - -(xi-ar 2 )a 3 /a; 1 +a4, d 5 = 
05, o?6 = <^7 = sxi. Hence D — 1- (— a%/x\) -a 3 - {[0:203 + £1(04 — a 3 )]/xi} -a;i-sx2 ■ sxi ~ 
— s 2 x 2 a^a^,[x 2 a^ + 0:1(04 — a 3 )]. The numerator is N = g 6 (2s) 3 , so the contribution from 5(b) 
is 



rp 9 6 s f If dxidx 2 

2ir z J a 3 as J x 2 [x 2 az + 2:1(04 — a 3 )J 

g- 6 slns /• ln(a 4 /a 3 ) 

= 2 / D q X i_F)qi^ ? r • 2.17 

2tt z J a 3 a 5 (a 4 — a 3 J 

The flow diagram 5(c) gives identical result so the total is 

T = T b + T c = -^V\/ 3 (A) , 
Ji(A) = [D qi±Dq2± (2.18) 

2.2.5 2nd, 4th, and 6th order QED diagrams 

QCD diagrams up to 0(g 6 ) are shown in Fig. 6, but for now we are interested only in those 
that contribute to QED electron scatterings. These are diagrams A, Bl, B2, C15-C20. 

Bl and C15 were computed above, others can all be computed in a similar manner. The 
contribution to T / (2s) from each diagram is || : 



.4 
Bl 

B2 

C15 



~9 2 h 
In (se~ 
2tT 



-9% 



+ ^9% 



2tt 
(Ins) 



2tt 2 



9 6 Js 
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C3 C4 C5 C6 C7 C8 




C9 CIO Cll C12 C13 C14 




C15 C16 C17 C18 C19 C20 




C21 



Figure 6: QCD scattering diagrams to 0(g 6 ). The QED diagrams come from the following 
subset: A, Bl, B2, and C15 to C20. 



C16 = -^VJ 3 



(Ins) 
~2^~ 



C17 = +^^ 6 (J3 + «/ 3 )=C18 



C19 = + ^V(J 3 - 7n/ 3 ) = O20 . (2.19) 

The total contribution up to 0(g 6 ) is obtained by adding up all the expressions in (2.19). 
The interesting, but sad thing, is that the (Ins) contributions all add up to cancel one 
another. It turns out that the 0(1) contribution from Bl + Bl is shown correctly in (2.19), 
but that of the sum of C15 to C20 is not. So in order to get a nonzero sum for those a much 
more difficult calculation accurate to 0(1) would have to be carried out. 



3 ABELIAN CUT DIAGRAMS 



3.1 Factorization Formula 

We want to deal with tree diagrams like Fig. 7, in which the momentum p is much larger than 
all the momenta q iy so that the following approximation for the denominator of propagators 
is valid: 

i i i 

j=l j=l j=l 

This is the same approximation used earlier in (2.9). 



15 



n 



qi q 2 q 3 q n 

Figure 7: A tree diagram, to be designated by the order of the dotted lines from left to right, 
as [12 • • • n]. 

A tree will be denoted according to the ordering of lines so the one shown in Fig. 7 is 
then [12 ■ ■ • n\. Its scalar tree amplitude is 



fn-\ 



a[12- 



n\ 



-2™y(5>i) MI 



3=1 



(3.2) 



In this subsection we will show that certain sums of amplitudes of this type can be expressed 
as a factorized product. 

We need some notations. If [Ti] are tree diagrams, then [T X T 2 • • • Ta] is the tree diagram 
obtained by merging these A trees together. For example, if [Ti] = [123] and [T 2 ] = [45], 
then [T]T 2 ] = [12345]. The notation {Ti; T 2 ; • ■ ■ ; Ta}, on the other hand, is used to denote 
the set of all tree diagrams obtained by interleaving the trees T 1; T 2 , ■ ■ • ,Ta in all possible 
ways. This set contains (J2a n aV-/ Y[ a n a- trees if n a is the number of gluon lines in the tree 
T a . In the example above, {Ti; T 2 } contains the following 5!/3!2! = 10 trees: [12345], [12435], 
[12453], [14235], [14253], [14523], [41235], [41253], [41523], and [45123]. 

Correspondingly, a{T±; T 2 ; • • • ; T^} will represent the sum of the amplitudes a[T] for every 
tree T in this set. The factorization formula Q states that this sum can be factorized into 
a product: 



a{Ti; T 2 ; ■ ■ • ; T A } = J] a[T a ] = a{T x \T 2 \ ■ ■ ■ \T A ] . (3.3) 

a=l 

Proof . Let [T] = [tit 2 ■ ■ - t n ] be a tree and a[T] its scalar amplitude defined by (3.2). Then 

a[T] 



To 



dn 



tn-l 



dTi 



tl 5 



(3.4) 



where the integration region T Q is defined by the ordering 00 > r tl > r t , 2 > • • • > r tn > 
—00. When summed over all T e {Ti; T 2 ; • • •}, the integration variables r ta retains only the 
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Figure 8: The abelian cut diagram [241 1 3 1 65] . 

ordering within each individual tree T», and for each tree they integrate from — oo to +oo. 
Using (3.4) again for individual trees T», we obtain (3.3). (end of proof ) 

There is a close affinity between the factorization formula and the string-like representa- 
tion. See Ref. [12] for a discussion on this point. 

The right-hand-side of (3.3) can be regarded as the amplitude a[Ti|T 2 | • • • \Ta] of the 
cut diagram [T x | T 2 1 • • - 1 T^] . This is just the diagram \T{T 2 - ■ - Ta] with cuts put on the 
propagators right after T 1; T 2 , ■ ■ • , T^-i- The cut amplitude a[Ti\T 2 \ ■ ■ ■ \Ta] is obtained from 
a\T{T 2 ■ ■ -Ta] by using the Cutkosky propagator —27rid(J2jUj) on the cut lines instead of 
the Feynman propagator (J2j + ^e) _1 - 

For example, the cut diagram [241 1 3 1 65] is shown in Fig. 8 and the corresponding cut 
amplitude is a[241|3|65] = a [241] a [3] a [65]. Note that the sections between cuts in a cut 
amplitude can be permuted at will, so for example, a[241|3|65] = a[3|65|241]. We shall refer 
to this as the commutative property of the abelian cut amplitude. 

The factorization formula can be used to sum up loop amplitudes as well. This is because 
(3.3) is valid whether are onshell or not, so the tree in Fig. 7 or 8 may very well be just a 
section of a much larger loop diagram, in which case (3.3) gives a sum of loop diagrams in 
terms of a single cut diagram. We shall call these abelian cut diagrams. 

Abelian cut diagrams resemble Cutkosky cut diagrams but they are quite different. Cuts 
here are only put on a high-momentum tree, and the resulting cut diagram represents a 
summation of Feynman diagrams rather than the discontinuity of one of them. 

In the simple case when every tree [Tj] = [i] is a vertex containing a single line, the 
factorization formula can be written as 

a[l; 2; • • ■ ; A] = a[l|2| • • - \A] . (3.5) 

This is called the eikonal formula and it has been known for a long time . 

A single abelian cut diagram represents a sum of Feynman diagrams, nevertheless it is 
easier to evaluate a single cut diagram than even a single Feynman diagram. This is mainly 
because there are fewer flow diagrams present in a cut diagram than a Feynman diagram. To 
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Figure 9: The six s-channel-ladder diagrams for n — 3. 



see that, suppose cuts are made on the tree of the top line, which is what we will usually do. 
A cut propagator is given by — 2iri5(sX) (see Sec. 2.2.2 for notation). Now y/sX is the ' — ' 
flow through the line, which must now be zero because of the presence of the 5-function, so 
a cut propagator is 'cut' also in the sense of severing the flow. This limits the various ways 
the ' — ' flow can go around the diagram and hence the number of possible flow diagrams. 



3.2 The Sum of QED s-Channel-Ladder Diagrams 

A (2n)th order s-channel-ladder diagram for electron-electron scattering is obtained by tying 
together the n photons from two electron trees like Fig. 7. There are all together n\ diagrams 
of this type depending on how the lines from the two trees are tied. Fig. 9 shows the 6 
diagrams for the case n — 3. 

These n\ diagrams can be labelled by the order of the photons \p\o-i • ■ • a n ] on the upper 
tree, when their order on the lower tree is fixed to be [12 • • -n]. 

Use (3.5) to sum over the n\ permuted upper trees, the result is an upper tree with all 
its propagators cut, as shown in the diagram between the two equal signs in Fig. 9 for the 
case n — 3. We may now use the commutative property of an abelian cut amplitude to 
symmetrize the photon lines attached to the lower tree. Use (3.5) again but now on the 
lower tree, we get a double-cut planar diagram, where all the electron propagators on the 
upper and on the lower trees are cut. A factor n\ is needed to convert the symmetrization 
into addition. For n = 3 this is shown in the last diagram of Fig. 9. 

The double-cut diagram is now easy to compute and one gets 



T 1 

Ys ~ ~n\ 



' n ~ 1 1 i 

]J(2s)g 2 Dq i± ——dqi + dq i ^(-2ni) 2 5(v r sq i+ )5(y/sqi 
l = \ Q>i o7T 



9 2 
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= -^-(-O^/nCA) , (3.6) 
n! 

where I n is defined in (2.13) and (2.15). 

Let us compare the result with those obtained in (2.19) by direct calculations. For n = 1, 
this should equal to A of (2.19) and indeed it does. For n = 2, it should be Bl + B2 and 
again it is. Note that individually Bl and B2 are of order Ins. To compute the sum which 
is of 0(1) correctly we need to know each Feynman diagram to sub leading- log accuracy. A 
cut diagram is already a sum so each can be computed just in the leading-log approximation 
as was done in (3.6). This ability to avoid subleading-log calculation is the main advantage 
of the cut diagrams. 

This is further illustrated by looking at n = 3, whose correct result is shown in (3.6) to 
be g 6 Iz/3\. This cannot even be obtained from (2.19) because the calculation of order g 6 
there is accurate only to leading log. 

In the impact-parameter space, (3.6) is given by 

— — = -i— -ig hip) . 3.7 
2s n\ 1 J 

Summing over all n, we get an impact-parameter representation of the T-matrix in the 
eikonal form: 

= -i J d 2 be tq ^{exp[-tg 2 h(b)] - 1} . (3.8) 

This elegant expression has been known for a long time ||. It is call the eikonal formula for 
high-energy scattering and is where the 'eikonal formula' (3.5) derived its name from. This 
S- matrix is unitary, with —g 2 Ii(b) being essentially the phase shift an electron accumulated 
when it shoots by another electron at an impact parameter b. High energy approximation 
enters by assuming the electrons travel essentially at a straight-line trajectory with fixed 
impact parameter (perpendicular distance between the scattered electrons) fixed. It leads 
to a total cross-section independent of energy. If pair productions in the intermediate states 
are included, then a rising total cross-section can be obtained 0. 



3.3 Other Sum Rules for QED 

In the last subsection we illustrated the power of the eikonal formula (3.5) in computing 
sums of all ladder diagrams. We shall give two examples in this subsection to illustrate the 
more general factorization formula (3.3) in computing other sums. These two examples are 
shown in Figs. 10 and 11. 

In each case, (3.3) is applied to the upper tree with the lower tree held fixed. The two 
cases correspond respectively to the identities 

o{31;2} = o[31|2] (3.9) 
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321 312 231 31,2 




123 123 123 123 



Figure 10: An example of using the factorization formula (3.3) to convert a sum of Feynman 
diagrams into a single abelian cut diagram. 




123 123 123 123 123 I 2 3 



Figure 11: Another example of using the factorization formula (3.3) to convert a sum of 
Feynman diagrams into a sum of abelian cut diagrams. 

a[321] +a[231] + a[132] = a{21; 3} - a{23; 1} + a{13; 2} 

= a[21|3] -a[23|l] +a[13|2] . (3.10) 

We shall now compute the abelian cut diagrams and compare them with the result 
obtained by computing Feynman diagrams. 

The flow diagrams for the four abelian cut diagrams are shown in Fig. 12. The flow in 
each case is unique, and the two poles for the two '+' integrations can all be taken at the 
bottom electron line. The computation is as before, and the result for Figs. 12 (a),(b),(c),(d) 
are respectively 

T a = - g %2s)(-i)±-(-lns)I 3 (A) = -T d 

T b = T c = 0. (3.11) 

The reason for the second equation is because in both 12(c) and 12(d) the arrows in one of 
the two loops are all in the same directions. 

The sum rules displayed in (3.9), (3.10), Fig. 10, and Fig. 11 demand that 

T 

C16 + C19 + C20 = (3.12) 
C16 + C18 + C17 = ±-(T b -T c + T d ) , (3.13) 

which can be seen to be true from (2.19), (3.12), and (3.13). 

In these two sum rules the leading In s is not cancelled, but the complicated transverse 
function J 3 in (2.19) is. 
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(c) 



(d) 



Figure 12: Flow diagrams to compute the right-hand-side of eq. (3.9) (diagram (a)) and 
eq. (3.10) (diagrams (b),(c),(d)). 



4 NONABELIAN CUT DIAGRAMS 



4.1 Multiple Commutator Formula 

In this section we consider the generalization of the eikonal formula (3.5) to nonabelian 
amplitudes. In other words, given a tree diagram (Fig. 7) with nonabelian vertices t a and 
amplitude a[12 • • • n] -t[12 ■ ■ ■ n] = a[12 • • • n] •t 1 t 2 • • • t n , we wish to find an expression for the 
sum of the n! permuted amplitudes to be expressed in terms of cut diagrams. This turns 
out to be the multiple commutator formula [|I[], which states that 

E «M*M = E ^Wo\ , (4.i) 

where S n = {1; 2; • • • ; n}, a [a] c is an abelian cut amplitude (see Sec. 3.1) for the cut diagram 
[cr] c , and t[a]' c is the corresponding nonabelian factor computed from the complementary cut 
diagram [a)' c . 

The cut diagram [a] c is obtained by putting cuts in the Feynman tree [a] in the following 
way. Proceeding from left to right, a cut is put after a number iff there is not a smaller 
number to its right. An external line is considered as a cut so there is never any need to 
have an explicit cut at the end of a tree. 

Here are some examples of where cuts should be put: [1234] c = [1|2|3|4], [3241] c = [3241], 
and [2134] c = [21|3|4]. 

The complementary cut diagram [a' c ] is one where lines cut in [a c ] are not cut in [o~' c ], and 
vice versa. Thus [1234]^ = [1234], [3241]' c = [3|2|4|1], and [2134]' c = [2|134]. 
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The nonabelian factor t[o~' c ] is obtained from t[o~] by replacing the product of nonabelian 
matrices separated by cuts with their commutators. For example, £[1234]' c = £[1234] = 
£ 1 £ 2 £ 3 £ 4 , £[3214]' c = £[3|2|4|1] = [£ 3 , [t 2 , [£ 4 , h]]], and £[2134]' c = £[2|134] = [£ 2 , £ x ]£ 3 £ 4 . 

Let me now illustrate this by writing down explicitly the multiple commutator formula 
for n = 3. In that case S3 = {1;2;3} contains 6 trees, and the right-hand side of (4.1) is 
given by 

a[l|2|3]£[123] + a[l|32]£[13|2] + a[21|3]£[2|13] + a[231]£[2|3|l] 

+a[31|2]£[3|12] + a[321]£[3|2|l] . ' "~ (4.2) 

As before, the trees in S n can be imbedded in much bigger Feynman diagrams, in which 
case the right-hand side of (4.1) and the corresponding nonabelian cut diagrams can be used 
to sum up loop diagrams as well. 

Application of this formula will be discussed in Sees. 5 and 6. Before doing that, we shall 
sketch the proof of (4.1) in the next subsection. 

4.2 Folding Formula 

The proof of the multiple commutator formula relies on another combinatorial formula which 
I call the folding formula [l[j : 

N N 

a[RoS] = ^2(-) k a{R;a 1>k .o}a[a k+ltN ) = ^2{-) k a{R;a ljk .o\a k+ i >N } ■ (4.3) 

k=0 k=0 

Let me explain the notation and the meaning of this formula. 

Consider the amplitude of a tree [RoS], put together from a tree [R], followed by a line 
[o], and then followed by another tree [S]. The tree [S] is assumed to have N lines, with 
[(Ti,fe] denoting its subtree formed from the first k lines, and a k +i,N its subtree formed from 
the remaining N — k lines. In other words, [S] = [o'i 1 fc fc+i,iv] f° r every k. The trees cr^o and 
&n+i,n are taken to be the null tree [0]. By definition, o[0] = 1. The notation cx ljfc means 
the tree o~i jk read in the reverse order. For example, if [S] = [14327856], then 0^3 = [143], 
a 4 , 8 = [27856], and cr 1>3 = [341]. 

Let [£] be a tree, with or without cuts. The symbol {T x ; T 2 ; • • • ; T A .t} is taken to mean 
the set of trees [Tt] for all [T] e {T x ; T 2 ; • • • ; T4}, and a{Ti, T 2 ; • • • ; T^.t} is the sum of these 
amplitudes a[T£]. 

The folding formula (4.3) shows how a[RoS] can be expressed as sums of products of 
amplitudes with line 'o' moved to the end of the tree in each case. Alternatively the line 
'o' is moved to a position just before a cut and by the commutative property of abelian cut 
amplitudes (see Sec. 3.1) this is just as good as moving it to the end. 

This formula is called the 'folding formula', or the cutting and folding formula, because 
the right-hand side can be obtained from the left-hand side first by cutting off the tree 
°~k+i,N at the end of [-RoS], and then folding the remaining tree about the point o, as shown 
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Figure 13: The folding formula illustrated. The cutoff part is not drawn. 

in Fig. 13. Finally, both branches of the folded tree should be interleaved to obtain the trees 
{R;a ltk .o} in (4.3). 

The folding formula is proven from the factorization formula (3.3) by an inductive and 
combinatorial argument I shall not repeat the somewhat involved arguments here, but 
I will work out some simple cases to illustrate the formula and its proof. 

Suppose we want to move line 1 of a [2 134] to the extreme right. Then in the notation 
of (4.3), we are faced with a situation where [R] = [2], [o] = [1], and [S] = [34]. Hence (4.3) 
gives 

a[2134] = a[21|34] - a{2; 3.1}a[4] + a{2;43.1} 

= a[21|34] - a[231|4] - a[321|4] + a [2431] + a[4231] + a[4321] . (4.4) 

Let me now show you how the proof goes for simple cases. Suppose [S] = [s] is a 
single line. Then (3.3) implies a[i?o|s] = a{Ro; s} = a[Ros] + a[R; s.o], which leads im- 
mediately to the folding formula a[Ros] = a[Ro\s] — a{R;s.o}. If [S] = [S1S2], then the 
proof goes in a similar way but is more involved. In that case we can use (3.3) to get 
a[i?o|siS2] = a{Ro; S1S2} = a[RosiS2] + a{R; S1.0S2} + a{R; S1S2.0}. The term a{R; S1.0S2} 
can be computed from the N = 1 formula, which yields a{R; Si.o}a[s 2 ] — a{R; Si, S2.0}. 
Substituting this back and noticing that 

a{R; S1S2.0} — a{R; si, S2.0} = — a{R; s 2 si-tf} , (4.5) 

we obtain the desired N = 2 result 

a[RosiS2] = a[Ro\siS2] — a{R; si.o}a[s2] + a{R; s^Si.o} . (4.6) 
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Similar proof goes for any N. 

I will now show how the multiple commutator formula (4.1) can be obtained from the 
folding formula (4.3). 

Take any abelian Feynman amplitude a[a}. Use (4.3) to cut and to fold, so as to move 
the number '1' to the end. The result is a sum of product of two a[- ■ •], the first one having 
the number '1' at the end, and the second one not containing the number '1'. Use (4.3) again 
to cut and to fold the second a[- ■ •], now to move the smallest number in its argument to 
the end, and so on. The final result is a complicated sum of products of several a[. . .], each 
having the smallest number within its argument located at the end. In other words, each is 
a cut amplitude a[p] c of some cut diagram [p] c determined in the way specified in Sec. 4.1. 

Conversely, given a [p] c , we can obtain all the [cr]'s that gave rise to it this way by repeated 
unfoldings and gluings. To see how this works in more detail, let us first consider the case 
when [p] c contains no cut. That means the number '1' is at the end, so it is of the form 
[p]c — b~l]- To unfold the tree [rl], pick out any two interleaving subsets [ti] and [f 2 ] from 
[t], in all possible ways, viz., find all such Tj so that [rl] 6 {n; r 2 .l}. Then the desired trees 
after unfolding are [a] = [71I.T2]. The sign involved in (4.3) is (— ) fc , where k is the number 
of boson lines in tree [r 2 ] . 

The nonabelian factors associated with a[a] is t[a], hence the nonabelian factor associated 
with a[p] c = a{rl] is XX - ) ^bi^], which is nothing but the multiple commutator t[p]' c . 

We have now considered the case when [p] c has no explicit cuts. If [p] c has explicit cuts, 
then we can apply the argument above to each of the cut sections, to obtain a multiple 
commutator for the nonabelian factor of each cut section. Thus the nonabelian factor for 
a[p] c is always t[a]' c as specified in Sec. 4.1, and the multiple commutator formula is proved. 
(end of proof) 

5 QUARK-QUARK SCATTERING TO 0(g 6 ) 

5.1 Color Factors 

The color factor of a quark-gluon vertex is t a , and that for a triple gluon vertex is if a bc- 
Putting these together, and using (1.8) and (1.9), or graphically Fig. 1, the color factor for 
a QCD Feynman diagram can be computed ||. In what follows we shall concentrate on 
'quark-quark' scattering where the 'quark' is allowed to carry any color in an SU(N C ) theory 
with any N c . In other words, the 'quarks' below could very well have been gluons as far as 
color is concerned. 

Fig. 14 contains an illustration as to how this can be done. For other examples see 
Ref. [6]. 

The color factor for a nonabelian cut diagram is computed in a similar way. Fig. 15 
contains one such example. 

To order 0(g 6 ), the color factors can all be reduced by this method to combinations of 
the four shown in Fig. 16 0. 
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= G 4 + 2 c G 2 - G 3 

Figure 14: An illustration for the computation of the color factor of a Feynman diagram 




Figure 15: An illustration for the computation of the color factor of a nonabelian cut diagram. 



G l G 2 G 3 G 4 

Figure 16: The four color factors for 0(g 6 ) diagrams. 
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5.2 Sum of Feynman Diagrams 



To 0(g 6 ) the QCD Feynman diagrams for quark-quark scattering are shown in Fig. 6. Each 
of their color factors can be computed by the method of the last subsection, and each of its 
spacetime amplitude can be computed as in Sec. 2.2. In fact, those diagrams common to 
QED have already been listed already in (2.19). For others, see Ref. [6] for details of the 
computations. 

Putting these together, The contribution from each diagram in Fig. 6 to T/(2s) is {(3 = 



A 


= -g 2 h-G 1 










Bl 


= -g 2 t3\n(se-™)l 2 
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= -g 2 (3 2 (\ns)2J 3 -(G 4 -G 3 + 


3cG 


2 + c'G,) 


C17 


= +g 2 (3 2 (\ns)(J 3 + mI 3 )-(G 4 


+ cG 


2) = C18 


C19 


= +g 2 (3 2 (\ns)(J 3 -mI 3 )-(G 4 


-G 3 


+ 2cG 2 ) = 



C20 (5.1) 
The functions J 3 and /„ were already given in (2.21) and (2.16). The function J 2 is 

d?q ± 1 



■/2(A) 



(5.2) 



The sum of these diagrams is 
1 



2s" 



-T 



-9 2 h 



+ 27 lnS 



1 — a In s + — a 2 \n 2 s 



■G 1 + l -(g 4 I 2 --g 6 I 3 \ns)-G 2 



■1 1 



■G 3 + ^ 6 /3-G 4 



where 



S(A) = ^A 2 / 2 (A) . 



(5.3) 



(5.4) 



There are several things worth noting: 
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1. In the fourth order, the leading term proportional to (Ins) is cancelled out between 
Bl and B2 in the color amplitude G 2 , though not in Gi. 



2. In the sixth order, the leading (Ins) contributions to G4 from C15 to C20 also add 
up to zero. The expressions given in (5.1) are not accurate enough to deal with the 
subleading terms as already noted in Sec. 3.2. The term in (5.3) proportional to G4 is 
obtained separately from the eikonal formula. 

3. As a result of these cancellations, the coupling-constant and energy dependences of the 
Gi amplitude is of the form g 2m (g 2 In s) p , where m — 1, 2, 2, 3 respectively for the color 
factors % = 1,2,3,4, and p is determined by the order of the perturbation. Referring 
to Fig. 16, we see that m is simply the number of vertical lines in the color factor, 
or as we shall see in the next section, the number of reggeons being exchanged. It is 
important to note that this dependence is not true for individual Feynman diagrams. 

4. To leading-log accuracy, the functions J 2 and J3 get cancelled out. Only I n appears in 
the final result. 

5. The quantity in the square bracket in the Gi term of (5.3) is simply the the first 
three terms in the expansion of exp(— a Ins). If we define the (transverse) reggeon 
propagator to be 

Rx (A, s) = h (A) exp(-a( A) In s) , (5.5) 
and analogous to I n = (*h) n we define the n-reggeon propagator to be 

R n (s,A) = (*R 1 ns,A) , (5.6) 
then the sum (5.3), to 0(g 6 ), is equal to 

= -g 2 Ri(A, s) • Gi + |^ [R 2 ( A, s) • G 2 + R 2jl (A, s) • G 3 ] + |^ 3 ( A, s) ■ G 4 . (5.7) 



The function i? 2j i is equal to 

1 



i* 2>1 (s,A) = V(lns) 

71 



-A 2 / 2 2 
2 2 . 



(5.8) 

The physical interpretation of these equations will be deferred till Sec. 6.1. 



5.3 Sum of Nonabelian Cut Diagrams 

The QCD nonabelian cut diagrams to 0(g 6 ) are shown in Fig. 17, with a subscript 'c' to 
emphasize that they are cut diagrams and not Feynman diagrams. Their computations can 
be found in Ref . [2] . 
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C16 c C17 c C18 c C19 c C20 c 

Figure 17: Nonabelian cut diagrams for QCD to order g 6 . 



For the s-channel-ladder cut diagrams, the computation is straight forward and one gets 
their contributions to T j (2s) to be 



Bl c = 

B2 C = 

C15 c = 

C16 C = 

C17 c = 

C19 c = 

C20 c = 



+ ±( ]nS )g%. cGl 

Ztt 
6 

" ; (lns)g 6 J 3 -c 2 G 1 



2tc' 
C18 c = 



__L,(l ns )^/ 3 .( c G 2 -G 3 ) 
Ztt 



(5.9) 



Similarly, we choose to cut the line of the planar diagram CI to obtain 

"1 



Cl c = -g 4 pi(lns) 
C2 C = -g 2 (3 2 (ln 2 s\ 



-A 2 /| - J 2 / 2 



A 2 / 2 2 - J 2 / 2 



G 3 
•c 2 G! 



(5.10) 



The twelve diagrams C3 to C20 can be divided into four groups of three, each giving 
identical contributions, so it is necessary to consider only one of these groups. The group of 
C3, C7, C12 have been chosen for that purpose. There is a symmetry between gluon lines 
2 and 3 so we may double this group and consider it as a sum of six Feynman diagrams. 
By applying the multiple commutator formula, the six cut diagrams shown in Fig. 17 are 
obtained. Their values are 



C3 r 



C7 r 



C12 r 



±(CX e + C3») = -g 4 (3i(\ns)^J 2 I 2 -G 3 
i(C7' c + C7' c ') = 

\{C12' C + C12" C ) = g 2 /3 2 (ln 2 s)^J 2 l2-{-c 2 Gi) 



(5.11) 
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It is easy to verify that they sum up once again to (5.3) as they should. 

Comparing this calculation of the cut diagrams with the earlier calculation of Feynman 
diagrams, one notes that 

1. The coupling-constant and energy dependence g 2m (g 2 In s) p , noted previously to be 
true for sums of Feynman diagrams under item '3' of the last subsection, is now true 
for individual nonabelian cut diagrams. In other words, delicate cancellations of the 
kind found in Feynman diagrams when they are summed never appear in cut diagrams. 
This is one of the main advantages of dealing with nonabelian cut diagrams. 

2. In leading-log approximation, functions J2 and J3 that appeared in individual Feynman 
diagrams but get cancelled in the sum never even occur here. 



6 MULTIPLE REGGEONS AND QCD 

6.1 The Reggeized Factorization Hypothesis 

Eq. (5.7) shows that QCD scattering takes on a simple form at high energies, at least up to 
0(g & ). The color factors present are those of Fig. 16, or equivalently Figs. 18(a), (b), (e), 
and (c). Interestingly enough, the spacetime amplitudes are also given by these same figures 
in Fig. 18, if they are interpreted as reggeon diagrams. 

The vertical lines in Fig. 18 are reggeized gluons. The horizontal lines are the lead- 
ing particles (top and bottom), or ordinary gluons produced from the reggeons. These are 
two-dimensional diagrams in the transverse-momentum space: the (vertical) reggeon prop- 
agators are given by Ri(s, q±) of (5.5), horizontal lines carry no factor, each vertex carries 
a factor g, and loop integrations for — T/(2s) are given by the measure —i J d 2 q±/ (2n) 2 . 
The spacetime amplitude for — T/(2s) from an m-reggeized-gluon (mrg) diagram without 
produced gluons (Figs. 18(a), (b), (c)) is then given by (— i) m ~ 1 g 2m (*Ri) m (s, A) /ml = 
(— i) m ~ l g 2m R m (s, A)/ml, where the factor 1/ml is a symmetry factor. This agrees with 
the terms Gi, G2, G4 in (5.7). 

For the reggeon diagram 18(e) with a produced gluon, the rule is more complicated. A 
Lipatov vertex gC^q^ q i+1 ) is to be placed at every reggeon-reggeon-gluon junction. For the 
elastic scattering diagrams in Fig. 18, they always come in a pair, so the factor is 



- 29 Qi+i) = 9 C^iqi, q i+1 )C^(A — q i: A — q i+l ) 



= -2g 2 



(A - q t± ) 2 q 2 +1± + (A - q l+1± ) 2 q, 



2 ' 

i± 



(qi± ~ qi+i±)' 



(6.1) 



as shown in Fig. 19. 
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(a) (b) (c) (d) 



+ 



+ 



(e) (f) (g) (h) 



+ 



+ 



(±) (j) (k) 



+ 



+ 



(1) (m) 



Figure 18: Reggeon diagrams for quark-quark scattering. 



1i 



A-q, 



A-q 



i+1 



Figure 19: (a) A Lipatov vertex gC '^(qi, qi+i); (b). Two Lipatov vertices —2g 2 IC(qi,qi 1 ). 
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In addition, the vertical reggeon lines should be thought of as covering the rapidity range 
stretched between the top leading particle and the bottom leading particle. Each horizontal 
line connected by two Lipatov vertices must be allowed to slide along this whole rapidity 
range to produce a factor (— 2ni) /2{2-n) 2 J dy ~ — Uns/in per produced gluon. In 0(g 6 ), 
the contribution from 18(e) to — T/(2s) is therefore given by 

2 

g 6 J d 2 q 1± d 2 q 2 ±Ri(s, q^Riis, q 2 ±)Ri(s, A - qi±)Ri(s, A — q 2± )- 



T 
2~s 



1 — i In s 



2! 4vr 
(-2)/C(g 1; g 2 
— ig 6 In s 



—i 



8vr 



—% 



A 2 



(2tt) 2 J J d2qi±d2q2± 
(A-q 1± ) 2 q 2 ± + (A 



1 1 



?L (A - <?i±) 2 (A - q 2 ±) 



(qi± ~ q2±f 



ig 6 In s 
2vr 



/, - -A 2 / 2 



(6.2) 



in agreement with the G3 term of (5.7). 

In the leading-log approximation where g 2 <C 1 and g 2 \ns = 0(1), every R4 and every 
/C are 0(1). Since the produced gluons must slide along the whole rapidity range mapped 
out by the reggeon lines, whenever one of them appears a factor g 2 \ias = 0(1) is produced. 
Hence an mrg amplitude is of order g 2m , whether there are produced gluons or not. This is 
certainly true in (5.7) but it should also be generally true for all reggeon diagrams in Fig. 18. 

This picture of reggeized exchange is simple and elegant, but does this 0(g 6 ) result 
generalize to all orders? In other words, can the sum of all high-energy elastic scattering 
Feynman diagrams be factorized into reggeon diagrams like those in Fig. 18? If it does, is 
the reggeon still given by (5.4) to all orders, or something more complicated is required? We 
shall refer to the affirmative answer to both of these questions as the reggeized factorization 
hypothesis. 

As discussed above, this hypothesis is completely verified to 0(g 6 ). It has also been 
partially checked to 0(g 8 ) and 0(g 10 ) M. The main focus of this section is to ask whether 
this is true to higher orders. To be sure, we do not know the full answer at the present 
but some partial answers are known. For example, using the technique of nonabelian cut 
diagrams, we can show that the hypothesis is correct at least for s-channel-ladder diagrams. 
This will be discussed in Sec. 6.4. Before doing that, we shall go over briefly what has been 
known for some time about multiple reggeon exchanges. 



6.2 BFKL Equation and Unitarity 

By summing t-channel-ladder diagrams in the leading-log approximation one obtains the 
result that the octet channel is dominated by a reggeized gluon exchange, with a regge 
trajectory given by a(t) = 1 + a(A), where a(A) is given in (5.4). The propagator of this 
reggeized gluon is given by (5.5). 0, [|. 
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The 2rg (two-reggeized-gluon) exchanges come from the sum of t-channel-reggeon-ladder 
diagrams || . These are diagrams in which n ordinary gluons are produced and absorbed by 
two reggeized gluons. They are illustrated in Fig. 18 by diagrams (b), (e), and (1) respectively 
for n = 0,1,2. 

In the leading-log approximation, the produced gluons are strongly ordered in rapidity, 
hence the rungs of the ladder do not cross. A diagram with n produced gluons has n + 1 
transverse integrations. The integrand consists of n factors of K, and 2(n + 1) reggeon propa- 
gators R\. It amplitude has been computed in the following way |§. First, its discontinuity is 
calculated from the Cutkosky rules. Then the whole amplitude is computed using s-channel 
dispersion relations. 

The color-octet amplitude for quark-quark scattering is predominantly real, and its lead- 
ing contribution is given by the lrg (one reggeized gluon exchange) diagram. However, this 
real part should also be obtainable from the imaginary part computed in the 2rg Cutkosky 
diagram via a dispersion relation. Explicit calculation shows that this is indeed the case. 

The color-singlet amplitude is the Pomeron amplitude originally proposed by Low and 
Nussinov flI5|. It is predominantly imaginary, a factor g 2 down from the color-octet am- 



plitude, and dominated by the 2rg diagrams in the leading-log approximation. Its com- 
putation is fairly difficult |5), and it leads to a total cross-section growing like S A , with 
A = g 2 (In 2) N c /it 2 . This violates the Froissart bound and unitarity, hence the computation 
as it stands must be incomplete. Subleading-log contributions must be added to unitarize 
it. These may be the subleading-log contributions from the same t-channel-reggeon-ladder 
diagrams, or the leading-log contributions from diagrams associated with more reggeon ex- 
changes, or both. 

Subleading-log contributions from Feynman diagrams are very difficult to compute, so if 
we must do it that way there would be very little hope of obtaining any result except in very 
low orders. However, if we sum up all multiple reggeon exchange diagrams as in Fig. 18, 
then formally s-channel unitarity is restored. From the discussion at the end of the last 
subsection, we know that an mrg amplitude is of order g 2m , so multiple reggeon diagrams 
computed to leading-log approximation do contribute to subleading-logs of the Pomeron 
amplitude. Thus there is a hope of unitarizing the Pomeron amplitude by including all these 
multiple reggeon diagrams computed to leading-log approximation, provided the reggeized 
factorization hypothesis is correct. 



6.3 Nonabelian Cut Diagrams vs Feynman Diagrams 

We saw in the 0(g 6 ) computation of Sec. 5 that delicate cancellations occur in the sum 
of Feynman diagrams, making it necessary to compute individual Feynman diagrams to 
sub leading accuracies in order to obtain a finite sum. In contrast, this never happens in 
nonabelian cut diagrams, so they can be computed simply to the leading order, a much 
easier task. We shall argue below that this situation is true to all orders provided the 
reggeized factorization hypothesis is correct. 

Recall that for nonabelian cut diagrams the color factor is computed from the comple- 
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mentary cut diagrams. For a diagram to contribute to an mrg amplitude, it must have m 
uncut propagators on the upper tree, which means that the corresponding cut diagrams from 
which the spacetime amplitude is computed must have m cut lines. 

Now an £-loop amplitude can grow at most like (Ins) . However, each cut propagator 
with its 5-function will deprive that loop of a In s factor because it replaces integrations like 
J £ dx/x ~ Ins by — 2ni J dxS(x) ~ 1. A diagram with m cuts is thus deprived of at least m 
factors of Ins, making it grow at best like g^2m)(g 2 In s) p , for a diagram of order 2(m + p). 
Since g 2 \ns = 0(1), this is of order g 2m , precisely what an mrg amplitude requires, so it is 
sufficient to compute each nonabelian cut diagrams just to the leading order. 

In contrast, the corresponding Feynman diagram has no cut in its spacetime part nor its 
color part. Its color factor will generally involve nrg diagrams with n > m. If this diagram 
contributes to an nrg diagram without delicate cancellation then its spacetime amplitude 
must be deprived of more (since n > m) factors of In s. On the other hand, without cuts 
in the spacetime diagram, there is no guarantee to have any deprivation of In s at all. The 
absence of cuts in both therefore conspire to render too many In s factors compared to those 
needed for reggeon amplitudes. In order for the reggeized factorization hypothesis to be 
valid, delicate cancellations are then necessary in the sum to reduce the In s powers. 

This argument shows that nonabelian cut diagrams are much more suitable for the com- 
putation of sums than Feynman diagram. In fact, there is virtually no hope of obtaining 
sums in higher orders by computing Feynman diagrams individually. 

6.4 s- Channel-Ladder Diagrams 

In this subsection we outline the computation of the s-channel ladder diagrams to show that 
they satisfy the reggeized factorization hypothesis. For details of the arguments see Ref. [10]. 
For reasons discussed in the last subsection, we will use nonabelian cut diagrams rather than 
Feynman diagrams in our computation. 

The Feynman diagrams are the same diagrams as those computed in Sec. 3.2 for QED 
to obtain an eikonal amplitude, leading to a constant total cross-section in that case. The 
situation in QCD is quite different. We know from the 0(g 6 ) computation in Sec. 5 that 
there are Ins factors present even for the Pomeron amplitude, so the total cross-section is 
no longer constant in QCD. 

There is another important difference. The sum of these ladder diagrams in QED is 
gauge invariant but that in QCD is not. To make it gauge invariant we virtually have to 
include all diagrams, and that computation is very difficult. What appears below is the 
result in the Feynman gauge. 

We shall consider quark-quark scattering in which the 'quarks' carry an arbitrary color in 
an SU (N c ) theory. Since spin is unimportant at high energies, what is thus obtained is valid 
for gluon-gluon scattering as well. The two quarks will be represented by two horizontal 
lines, and the n exchanged gluons by vertical or slanted lines as in Fig. 20 for n — 4. 

We shall choose the planar diagram to be the one in which all propagators on the upper 
line are cut, as shown in Fig. 20(a). This can be accomplished by numbering the gluons in 
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the manner shown. For convenience we shall also draw the gluons of the planar diagram as 
vertical lines. 

All the other nonabelian cut diagrams can be obtained by pulling the upper ends of a 
number of gluon lines leftward, in all possible combinations. The propagator to the right of 
a slanted line is uncut, and the propagators to the right of a vertical line is cut. See Fig. 20. 
This is so because once the upper end of a gluon line is moved leftward, there will always be 
another gluon with a smaller number to its right. 

We shall refer to these s-channel-ladder cut diagrams as SC diagrams. Their comple- 
mentary diagrams used to compute color factors will be referred to as SCC diagrams. 

As discussed in the last subsection, a nonabelian cut diagram with m cuts are deprived 
of at least m would-be In s factors. Such a diagram is said to be saturated if it is deprived 
of exactly m factors of Ins, and no more. The leading contribution to a mrg diagram will 
come only from the saturated SC diagrams with m cuts. 

By referring to the computations in Sees. 2.2.3 and 2.2.4, we see that even diagrams with 
no cuts may not get its full share of Ins factors - Fig. 4 does but Fig. 5 does not. Like 
the cases in Figs. 4 and 5, an SC diagram is saturated iff it does not have two consecutive 
uncut propagators on the upper line. Thus 20(k) is not saturated, but every other diagram 
in Fig. 20 is. 

We outline now the procedure for computing sums of SC diagrams of any n. For details 
please consult Ref. [10]. 
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We shall refer to color factors of the type shown in Fig. 18 as regge color factors. Given 
an SCC diagram, the first step is to use eqs. (1.8) and (1.9), or equivalently the graphical 
rules in Fig. 1, to decompose its color factor into a sum of regge color factors. We will 
call a regge color factor primitive if it remains connected after the top and bottom lines are 
removed. Thus the regge color factors in Figs. 18(a), 18(e), 18 (k), and 18(1) are primitive, 
but none of the others explicitly shown in Fig. 20 are. 18(1) and 18 (m) do not occur in SCC 
diagrams so we will ignore them. The other three primitive ones 18(a), (e), and (k) will 
respectively be named Hq, Hi, and H 2 . It can be shown that the only primitive regge color 
factors entering into SCC diagrams are H p (p = 0, 1, 2, 3, • ■ ■), representing a pattern with p 
horizontal steps climbing up from left to right. 

It can also be shown that regge color factors with the same primitive components can be 
identified in the leading-log approximation. For that reason a regge color factor can simply 
be labelled by the number of each primitive component it contains, i.e., can be expressed in 
the form $ = n^Lo ■ For example, the color factors of 18(b), (c), (d) are Hq, Hq, and Hq 
respectively. The others are 18(/) = 18(g) = H H U 18(h) = H$H U I8(i) = 18(j) = H\. 

The next step is to collect all SC diagrams contributing to a given regge color factor $. 
It turns out that these diagrams can be summed up using the factorization theorem (3.3) to 
be (A = -T/2s) 

^{$}(A) = nH) /p_1 (A) , (6.3) 

where A{H p } is the sum of all SC amplitudes whose corresponding SCC diagrams contain the 
connected primitive color factor H p . This shows that the reggeized factorization hypothesis 
holds for all s-channel-ladder diagrams. 

In particular, A{H }(A) is simply g 2 R\(s, A), where R\ is given by (5.5) to order g 2 . The 
g A contribution comes from diagrams Cl c and C2 C of Fig. 17 and they are not SC diagrams. 

Using (5.6), the G2 and G4 amplitudes of (5.7) are special cases of (6.3) with f = 2, 3 
and all other f p — 0. 

Similarly, A{Hi}(A) = (—i/2)R 2t i(s, A), where R 2t \ is given by the first term of (5.8) - 
the second term there does not come from SC diagrams. 
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